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Abstract 

With a modest revision of the Standard Model, in the fermion mass sector only, the systematics of 
the fermion masses and mixings can be fully described and interpreted as providing information on 
matrix elements of physics beyond the Standard Model. A by-product is a reduction of the largest 
Higgs Yukawa fine structure constant by an order of magnitude. The extension to leptons provides 
for insight on the difference between quark mixing and lepton mixing as evidenced in neutrino 
oscillations. The substantial difference between the scale for up-quark and down-quark masses 
is not addressed. In this approach, improved detail and accuracy of the elements of the current 
mixing matrices can extend our knowledge and understanding of physics beyond the Standard 
Model. 
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I. INTRODUCTION 


For several decades, the quantum numbers and corresponding gauge interactions that dis¬ 
tinguish the different generations of fermions have been sought without overt success. The 
various efforts to understand the fermion masses have ranged from substructure to Grand 
Unihcation. The former approaches include rishons[T] and technicolor [2], but always en¬ 
counter a fundamental problem: For relativistic constituents, the limiting gaps between 
eigenstates tend to a constant value, as in the MIT bag model|3] (surprisingly parallel to 
the eigen structure for the non-relativistic harmonic oscillator). With sufficiently strange 
potentials (e.g., the “dracula” potential |1]), the lowest few states may be forced to match 
the increasing gaps found in the real world, but they ultimately tend to a constant gap and 
predict additional states within the range of experiment that remain unseen. The latter 
approach has found some interesting relations between quarks and leptons O |6] and even 
some mixing angles [7] but no convincing overall solution has been obtained. Efforts along 
these lines have continued. [Sj 

We consider a revision of the Standard Model (SM), in the fermion mass sector only, by 
taking a contrary point of view, namely that, within the (revised) SM, all of the fermions 
with a given electric charge should be viewed as having nothing that makes their right- 
chiral, weak interaction singlet components distinguishable to the Higgs boson. That is, we 
take the Higgs coupling to be completely insensitive to “generation” and discard the Yukawa 
coefficients that have been (artihcially) inserted in the SM to reproduce the observed masses. 

In the past, many others have constructed similar mass matrices by starting from various 
symmetry assumptions. A few of them, of which we are aware, are referenced here. imni 
They are all “top-down” approaches, making initial symmetry assumptions. Conversely, we 
take a “bottom-up” starting point focusing on the absence of a known symmetry or quantum 
numbers. This may be viewed as an accidental S 3 symmetry but that is not the fundamental 
nature of our basic assumption. Rather, we assume that the hnal form of the fermion mass 
matrix is determined by all possible (loop) corrections to the SM from physics beyond the 
SM (BSM physics) and invert the relation to extract information on some matrix elements of 
BSM physics. This should be contrasted especially with the approach in previous work, such 
as for example in Refs. m), where the number of unknown BSM parameters is reduced by 
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assuming only limited forms for the initial mass matrices. 

We apply this concept here to the quarks and comment on related implications for leptons, 
reserving a full discussion of the leptonic system, including neutrinos and their additional 
complications, to a later paper. This approach amounts to a small change in the mass 
sector of the SM, (hence, a modest revision, albeit a “massive” one) which provides an 
interpretation, in terms of BSM physics, of the deviations from the direct mass implications 
of our revised version of the SM. 

We obtain a consistent description of the quark masses and weak interaction current mix¬ 
ing for BSM corrections on the order of 1%, and constrain relations among the (many) 
allowed BSM parameters, with values that are “natural” in the classic sense (0(1)). We 
hnd that non-Cartan sub-algebra components of the BSM symmetry-breaking corrections 
are required. 

Although this approach proposes a resolution of the differences in fermion masses between 
“generations”, it does not attempt to address the differences in mass scale between the 
fermions of different electric charge, generally described as within “families”. Our approach 
does, however, reduce the number of these questions to two for the quarks and provides for 
a similar reduction for the leptons. 


A. Mass in the SM 


We hrst briefly review the structure of SM mass terms that we will revise. They are set by 
(arbitrary) Yukawa couplings of the (one, or in supersymmetry, at least two) Higgs boson(s) 
to the weak interaction active left-chiral doublets and weak interaction “sterile” (except for 
their U{1)b quantum numbers) singlets. The form is 


= E,;Yr„: ( Tf/, < 0 °, 0 + > ^ ^ 




'^Di 

d'Di 


-|- h.c. 


+ h.c. 


( 1 ) 


where the Higgs may be the same or different in the two sets of terms. 
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Unlike the comparable terms for the interaction with the weak gange bosons, which is natn- 
rally in the cnrrent basis, this Dirac notation is in the mass eigenstate basis and suppresses 
the information that pairs of independent Weyl spinors are involved. In the SM, the Yi are 
taken to equal the mass of the appropriate fermion divided by the vacuum expectation value 
of the Higgs boson. 

In a more general notation, a Dirac bispinor is composed of two Weyl spinors m 

T ^ 

qj = 

X 

where these two chiral spinors are constructed from another pair {(a and (b) with a hxed 
phase relation, as 


( 2 ) 


^ — Ca + Cb ( 3 ) 

X = ^2{-Ca+Cb) ( 4 ) 


SO that the Dirac mass term appears as 

^4^ = -(x^e + e^x)- 


( 5 ) 


This makes it clear that the Dirac mass term couples two independent Weyl spinors (left- 
and right-chiral) one of which is weak interaction active and the other sterile as described 
above. 


For completeness, we display the Majorana form for a Weyl spinor and its mass term: 


Tm = 






( 6 ) 


which will be relevant when we refer to neutrinos later. 


B. Basic concept 

The weak iso-singlet produced by the coupling of the Higgs boson to the left-chiral projection 
of the quark (generally, fermion) doublet has no known quantum number besides that of the 
U(l) that completes the electric charges of the helds. Similarly, the right-chiral projections 
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have no known other quantum number. Hence we are free to choose bases in the left- and 
right-chiral representation spaces such that the entries to the mass matrix for all of the quark 
(fermion) helds of a given charge are identical. The corresponding mass matrix appears to 
have been hrst described as “democratic” by Jarlskog [lO], although, for larger dimensional 
forms, it is known in nuclear and condensed matter physics as the origin of the “pairing gap”. 
The eigenstates for such a system consist of one (non-zero energy or) massive state with all 
other (two in the case of interest here) eigenstates being (at zero energy or) massless. This 
provides for a natural starting point consistent with the large gap between the heaviest and 
the lighter known fundamental fermions (of each charged type). We address the question of 
basis choice below in SeclIIIl 

Experimentally, however, the lighter fermions in each grouping are not massless. This de¬ 
viation from zero mass for the lighter fermions can be accommodated by assuming small 
deviations from “democracy”, consistent with perturbative corrections from BSM physics, 
which quite naturally allows for small mass eigenstates in the hnal result. A particular 
constraint on the parameters of the BSM physics is required for one of these to be much 
smaller than the other. Under our assumption, these effects afford a glimpse into the nature 
and structure of the BSM physics by requiring specihc relationships between some matrix 
elements (that would be calculable from any given BSM theory). 

The constraints can be discerned by examining the matrix, U, that relates the BSM-corrected 
mass eigenstates to the weak interaction eigenstates, for the up-quarks, and the correspond¬ 
ing matrix, V, for the down-quarks, combined into the form that produces the Cabibbo- 
Kobayashi-Maskawa [131 E] (CKM) matrix in the form as described by the Particle Data 
Group [I5] (PDG). That is, our approach recognizes that the Higgs coupling to the “active”, 
left-chiral, weak iso-doublet fermions is aligned with the weak gauge boson coupling, and so 
is diagonal in that (Weyl spinor) basis. As stated above, however, we take as a given that 
the right-chiral, weak iso-singlet Weyl spinors all present themselves indistinguishably to 
the weak iso-singlet object formed by combining the Higgs doublet with the active fermion 
doublet. The GKM matrix is produced conventionally by the misalignment between the 
mass and weak eigenstates, but here determined from the two separate components, U and 

ut. 

Another feature of the mass eigenstates of the 3x3 “democratic” matrix is that the eigen- 
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vectors for the mass eigenstates can take the form of a tri-bi-maximal (TBM) mixtnre of 
the original (cnrrent) eigenstates. One state mnst be maximally mixed; for the other two, 
one has a degeneracy choice. However, there is no impediment to an overall TBM choice. 
While nsefnl to simplify calcnlations, this is not particularly significant for the quarks, as 
the TBM structure cancels out in construction of the CKM matrix. (See Eqs.( 49p0 ) be¬ 
low.) However, it does have signihcant implications for the difference between the CKM 
and the corresponding Pontecorvo-Maki-Nakagawa-Sakata |T6] (PMNS) matrix for mixing 
in the lepton sector. We will comment upon this difference in our conclusions. 


II. STARTING POINT AND POSITIVE INDICATIONS 


The tri-bi-maximal (TBM) matrix 



1 

1 

1 


P6 

P2 

V3 

TBM = 

1 

1 

1 

P6 

V2 

V3 


2 

0 

1 


. V6 

^/3 . 


diagonalizes the “democratic” matrix 


Mdem — 2 X 


111 

111 

111 


to 


Mm = TBM^ X Mdem X TBM 
0 0 0 
= 0 0 0 
0 0 1 


(7) 


( 8 ) 


(9) 


where we have chosen the overall scale so the nonzero eigenvalue is unity. 

The efficacy of the “democratic” conjecture can be tested by inverting the TBM transforma¬ 
tion on the known quark masses (taken from the PDG|T5] and) placed into diagonal mass 
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matrices, viz., 


0 


and 


ruu 


2.3 0 

0 1275 0 

0 0 173500 


rrid 


3.8 0 0 

0 95 0 

0 0 4150 


( 10 ) 


( 11 ) 


where all values are expressed in MeV/c^. (We will ignore the signihcant uncertainties and 
variation with scale of these masses [T7] as the ratios vary less dramatically, although the 
values of even the ratios are not known to very high accuracy.) 


We now transform these inversely using the TBM matrix given in Eq.([^ and find that the 
resulting mass matrices are indeed quite “democratic”: 


l^u-TBM 


TBM X X TBM^ 


(173500) X 


0.33701 0.32966 0.33333 
0.32966 0.33701 0.33333 
0.33333 0.33333 0.33334 


( 12 ) 


and similarly 


rrid-TBM 


(4150) X 


0.34493 

0.32204 

0.33303 


0.32204 

0.34493 

0.33303 


0.33303 

0.33303 

0.33394 


( 13 ) 


where we have scaled out the overall factor of the largest mass in each case. Although the 
true accuracy is, of course, far less, we keep the extra digits to display which matrix elements 
are not identical after the (inverse) TBM transformation and so convey the patterns that 
will survive even substantial (within experimental uncertainties) changes in the ratios of the 
diagonal values. 


This demonstrates that only perturbatively small corrections to a democratic starting point 
are needed. (We have ignored CP-violation considerations here, but will return to them 
below.) The deviations from “democracy” are exceptionally small, less than 1% in the up- 
quark sector and less than 4% in the down-quark sector (for positive and negative deviations 
from an average). It is clear from this that something close in structure to Mdem (times an 
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overall mass scale, m) is a reasonable ansatz to consider for an initial mass matrix. (A 
similar result holds for the charged leptons.) 

This result conhrms that the wide range of quark masses is well described by an almost 
“democratic” mass matrix for each charge set of quarks, leaving only the overall scale differ¬ 
ence between up-quarks and down-quarks (and also leptons) to be understood. We do not 
address that difference here. 


III. A MODESTLY REVISED MASS SECTOR 


With these comments and results in mind, we propose the mrSM (modestly revised Standard 
Model) which differs from the SM only in the fermion mass sector. In terms parallel to those 
of the SM used above, we have 


^mr 






1 + 75 


^Di < > 


1-75 

2 

1-75 


_ 

Tdj 


-f h.c. 


-I- h.c. (14) 


There is now only one Y each for all of the up- and down-quarks and these two values are 
approximately one third of the value of the two largest in the SM. 


Of course, as long as each up-quark is related to a particular down-quark by the weak inter¬ 
action, one may choose a “rotation” of the pair to any basis. We reiterate the basic points: 
Because there are no known quantum number restraints, we are free to rotate the bases for 
both the left- and right-chiral helds independently to ensure that the mass matrix is demo¬ 
cratic. (With some particular BSM theory, this should be the basis required by quantum 
number constraints.) In that basis, it is clear that BSM corrections are perturbatively small, 
a situation in physics always much to be desired. This is especially so here, since effects of 
BSM physics, if they exist at all, must be suppressed. 


We note, as lagniappe, that this factor of 3 reduction in these remaining two Yukawa cou¬ 
plings to the Higgs held from the largest value in each case in the SM improves the formally 
perturbative character of this part of the SM by reducing the (set of) Higgs hue structure 
constant(s), Y^/dvr, by an order of magnitude. Of course, the total size of the effects of 



quark loops is unaltered, as the sum over all channels produces the same signihcant net 
effect as in the diagonal mass basis. However, the size of the total contribution becomes due 
to the number of diagrams contributing, not to any individual large one. 

We take the full Higgs plus BSM-loop-corrected mass matrix to have the form 


M-mrSM = m X [Mdem + ^M-BSm] 


( 16 ) 


now in the current quark basis consistently defined by the Higgs and weak vector boson 
couplings. That is, we define the mass matrix for each set of 3 quarks of a given electric 
charge as m (an overall scale which is approximately one-third of the mass of the most 
massive of each triple of the fermions of a given non-zero electric charge) times the matrix 
A4, where 


M 


(1 + eO 


X 


111 

111 

111 


\/lyo + 2/3 + :^y8 2/1 - ^ 2/2 2/4 - ^ 2/5 

2/1 + Iy2 \/lyo - 2/3 + ^ 2/8 2/6 - lyj 

Va + Iy5 2/6 + lyi V^2/o - ^2/8 


(16) 


This allows for the most general set of perturbative (for small e) deviations possible for 
a Hermitean 3x3 matrix from the democratic mass matrix produced by uniform Higgs’ 
coupling in each quark charge sector. The coefficients are chosen to match the normalization 
of the standard Gell-Mann SU{3) {U{3)) basis matrices. 


The BSM corrections are all taken to be proportional to the small quantity, e, defined by 
the diagonal matrix of known mass eigenvalues, (again, with the overall scale factored out) 


m X 


6(5 0 0 
0 6 0 
0 0 1 


( 17 ) 


The overall factor (1 -|- e.^) in Eq.(16) is introduced to rescale the largest eigenvalue of Ai to 
unity, needed to account for the effect of the eyo correction from the BSM physics. We will 
see below how a correlation between some of these Hi constants (which we view as matrix 
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elements of BSM physics loop corrections, see below) reduces the smallest eigenvalue from 
e to e6. 

As is apparent from Eq.([l7|, 5 is the ratio of the lightest mass of the three quarks (with the 
same electric charge) to the mass of the intermediate mass quark, and e is the ratio of that 
quark to the most massive of the three. In particular, for the quark mass values referred to 
above. 


e„ = 7.35 X 10■^ 5u = 1.8 x 10"^ 
ed = 2.29 X 10"^ 5d = 4.0 x 10"^ 


(18) 


Even the largest of these values easily qualihes as a small expansion parameter. We will see 
below that the 5s do not signihcantly influence our results, so the largest perturbation is 
provided by e^. 

The quantities yo, and ys describe a subset of the possibilities for symmetry breaking from 
the BSM physics: While (1 + e^) is only an overall scale revision, y^ and ys parameterize the 
usual Cartan sub-algebra of SU{3) (dynamical) symmetry breaking allowed for 3 complex 
degrees of freedom, and yo parameterizes the correction allowed by the U{1) factor of the 
overall U{3). These are all multiplied by a factor of e in the expectation that the BSM 
corrections are small, as is suggested by the data referred to above. It is apparent from the 
numerics above that this assumption is self-consistent. 

We also allow for off-diagonal corrections, reflecting a possible more complete breaking 
of the (accidental) U{3) symmetry. We have included these elements of necessity, having 
found that, with only ys and ys, the desired result of the CKM matrix (in the form of 
the combination UV^) cannot match all of the (moduli of the) elements of that matrix as 
reported by the PDG [IS]. However, as emphasized by Leviatan [TS|, some elements of 
a partially broken symmetry may survive when additional components outside the Cartan 
subalgebra, or even outside the parent spectrum generating algebra, appear. Interestingly, 
the converse does not hold. We will see below that it is possible to describe the mass and 
current misalignment with at least one and possibly both ys and ys vanishing as long as 
components outside the Cartan subalgebra contribute. 
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IV. EXPLICIT CONSTRUCTION 


After TBM transformation as in Eq.([^, an intermediate form of the mass matrix with BSM 


corrections as given in Eq. (16) becomes 


Ad ini 


0 0 0 

0 0 0 

0 0 1 + e^ 


+ e X 


( 19 ) 


Ivo 

^ (2/3 


\/3 
2/8 + 


§2/8 + |(2/i - 22/4p6) -^(2/: 

VirnG -^(2/2 y^ml)^ 

;^(2|/l — 2/4p6) + ~^Iybp7 ~ 


' (22/3 + 


■+6 


yim% + -^( 2/2 — ybmr)) 
^2/8 - 2/1 
2/4m6 + I (22/2 + 2/5m7) 


§ 2/8 + 


3 - 2/0 + 


3-- ■ 3V2^‘^y^ 

;^( 22/3 + yAm& I {‘2y2 + ybrn?)) 

§2/0 + §( 2/1 + 2/4p6) 


We have introdnced a p and m simplifying notation: 1 / 4 ^ — i/q = ?/ 4 m 6 , 2/4 + 2/6 = 2 / 4 p 6 , 
y 5 — yj = y^mi and y^ + yj = 2 / 5 p 7 to indicate snms and differences of the nnmerically 
labelled i/j’s. This combines them into shorter labelled forms to reduce the complexity of 


the appearance of the form of Aiint- We now reformnlate Eq.(19) in terms of phases as 


M, 


§ 2/0 


0 0 0 

0 0 0 +e X 
0 0 1 

- + 3^8 + 1 ( 2/1 - 22/4; 
-yae'"^ 
ybc-^y 


ipd) yaf' 

\y^ + ^2/8 ~ 2/1 


2/be 


1/3 


z,-n 


-yce 


i-y 


e + 


-yce 
2/0 + §( 2/1 + 2/4p6) 


( 20 ) 


where 


ya y^md)'^ + (2/2 2/5m7)^ 


2/2 — 2/5m7 
,2/3 — 2/4m6, 


a = — arctan 
2/6 = (2\/3i/ 8 + 2l/i - l/4p6)^ + 9(l/5p7)^ 


/3 = arctan 


32/5p7 


,2V3!/s + 2!/i - Pipe, 


( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 
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( 25 ) 


Uc — ^^\/( 22/3 + VAmQY + ( 2?/2 + UbmlY 
, ( 22/2 + 2/5m7 \ 

7 = arctan -- 

\ 2j/3 ~\~ / 


(26) 


and we take the positive values of the square roots. 

There is a phase freedom associated with each of the three fermions or eigenvectors. Al¬ 
though the overall phase can have no effect, we can make use of two phases to transform 
the mass matrix to the form 


Me 


0 0 0 
0 0 0 
0 0 1 


+ e X 


§ 2/0 - ^ 2/8 + 1 ( 2/1 - 22 / 4 p 6 ) 


2/a 


-2/a 




-2/c ^ + 




Vb e 

-2/c 

2/0 + liVi + Vips) 


(27) 


where 

C = a + /3 + 7 (28) 

is the only physically meaningful quantity and 2/a, Vb and 2/c are unaltered. 

We have come this far by transforming the mass matrix by the TBM matrix as in Eq.(|^, 
but we need further transformations to carry out the diagonalization. Here, we only display 
this calculation through order e. We proceed in two steps, hrst block diagonalizing with a 
generic matrix, A 3 _,. 2 , which will apply for both up-quarks and down-quarks with appropriate 
parameter values. We then complete the calculation by diagonalizing the remaining 2x2 
block with the generic matrix, X 2 a; 2 - The total transformation is given by 


Xtot — TBM X ^ 3^2 X ^2x2 


(29) 


We use the notation, X, to reflect the fact that the form of the matrices that produce the 
diagonalization of the (approximate) mass matrix will be applied to both the up-quarks, 
where the full matrix that produces diagonalization is conventionally labelled as U, and to 
the down-quarks, conventionally labelled as V. 
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A. Block diagonalization 


We now first block-diagonalize the 3x3 matrix in Eq.(27). With the expectation that the 


third component will dominate the eigenvector for the large (unit) eigenvalue of the full 
system, we choose an eigenvector of the form 

a 

vec3 = 


1 


(30) 


where we expect a and (3 to be of order e, and solve the hrst two of the equations in 

vec3 = 1 • vec3 (31) 

for a and (3. Their values to leading order in e are 


a = ey be 

13 = -eye 


< 


(32) 

(33) 


The value of ^ is determined by the requirement that the trace of equal the trace of the 


scaled matrix in Eq.(17): 


f = 1 + 3 - V^yo 


(34) 


Then the value of y^ can be found from the third component relation in Eq.(31). To the 

leading order approximation in e used above, it is 

1 


2/0 — + (5) + 2(?/i + |/4p6)) 


(36) 


It is now straightforward to construct two vectors orthogonal to vec3. With these, and after 
normalizing all three vectors, we can block diagonalize as an intermediate step to the 
full diagonalization. Since we are only working to order e, we do this in a way that minimally 
affects the 2x2 subspace, choosing 


vecl 


1 

0 


—a 


vec2 


0 

1 


-p 


(36) 


(37) 
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The minus signs appear because the full matrix of eigenvectors is necessarily of the form of 
a complex rotation and (3 is real, so it does not require complex conjugation. 

Together with uec3, the two vectors vecl and uec2, allow us to present an explicit unitary 
transformation matrix, X 3 _,. 2 , that produces, through order e, a block-diagonalized form of 
namely 

1 0 e i/fe 6 *“^ 

A3_s.2 = 0 1 —euc (38) 

_-eybe~^^ eye 1 _ 

which block diagonalizes to the desired unit eigenvalue and a 2 -by -2 block 

i “ 73^8 + I ~ -Va 

-Va |(1 + 5) + ^ys - \{2yi - y^^po) 

(The first two entries of the third row and column of the full matrix are zero to order 
by construction; the values of ^ and yo above may be used to conhrm the unit value of the 
resulting (3,3) element.) 

Note that, due to the isolation of the phase to the (1, 3) element of there are no order 
e imaginary terms in M 2 x 2 - It is here that the advantage of avoiding an order e imaginary 
contribution in M 2 x 2 becomes apparent: Once the factor of e is removed from M 2 x 2 , it is 
clear that X 2 a; 2 , the matrix that diagonalizes M 2 x 2 i would include an eigenvector with an 
order one imaginary phase which would lead to a large CP violation, inconsistent with 
observation, and so require a delicate cancellation to occur. The isolation of the net phase 
to the (1, 3) element of enforces this necessary result automatically. (A parallel result 
obtains if the phase freedom is used to combine the separate phases in Al^ into the invariant 
total phase (, but in the (2,3) element of Al^ instead.) 




B. Completion of Diagonalization 


To complete the analysis and proceed to apply constraints on the BSM parameters, we 
need to diagonalize M 2 x 2 - This is, of course, straightforward; we simplify the notation by 


identifying 

1 1 

Vf = - 3 ( 21/1 - y^pe) 


(40) 
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so that the eigenvalues are 


Next we dehne u by 


and set 


^(1 + ^) ± \Jyl + y} 

tan(2a;) = 

yj 


ya 

yf 


1-5 

2 

1-5 

2 


sin(2ci;) 


cos(2a;) 


(41) 

(42) 


(43) 

(44) 


SO that the required values of the eigenvalues, 5e and e, are obtained. The 3x3 matrix that 
diagonalizes the 2x2 block in the full mass matrix is then just 




cos(a;) sin(c(;) 0 
—sin(a;) cos(a;) 0 
0 0 1 


(45) 


where we have relabeled X 2 x 2 as for notational simplicity. The full matrix that completely 
diagonalizes the mass matrix to order is now 


Xtot = TBM X X3^2 X X, 


(46) 


as we described in Eq.(29). Explicitly, 

cos (a;) 


Xtot = 


-sm a; 


sin^caj 

cos(a;) 


m e 


< 


-ei/c 


-€{llb€^cos{uj) + ycsin(w)) -e(y,e'^siii(w) - y^cos(w)) 


(47) 


While every particular BSM physics model will specify the values of the that enter into the 
determination of u, it should be clear from the above that inversion of the mass spectrum 
is not sufficient to determine a value for u. Thus, it remains an effective free parameter in 
what follows. 


V. FITTING TO THE CKM MATRIX 


To compute the CKM matrix, we need the result in Eq.(47) evaluated for both the up 


quarks, and for the down-quarks. Again, the separate matrices for these are conventionally 
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labelled U and V respectively [15], so that 


CKM = U xVl (48) 

However, we reiterate that the PDG description is one in which these matrices transform 
from mass eigenstates to current eigenstates, but our derivation above is for the transfor¬ 
mation of current eigenstates to mass eigenstates. Hence, the Hermitian conjugates are 
interchanged and the of the PDG is our Xtot for the down-quarks and similarly, their U 
is the Hermitian conjugate, of our Xtot for the up-quarks. 

On combining the results for up-quarks and down-quarks to produce the equivalent of the 
GKM matrix, 

CKM = XP,. TBAft TBM X^_ (49) 

we see that the TBM factor cancels out in the product. Hence, it is sufficient to calculate 

CKM = A'i. X A'P,, X X3^2, X (50) 

Gontinuing to work only at hrst order in the small quantities, we dehne 

Adu = edVbde^^'^ - euVbue^^^'^ (51) 

Bdu i^^dUcd ^uUcu) 

and again make use of phase freedoms: In the GKM matrix, there are six available, three 
each from the up-quark and down-quark sectors. As before, one overall phase can have no 
effect, but four of the remaining hve can be chosen to obtain 

GAM(1,3) = [cos(c<;„)Ad„ - sm{uJu)Bdu\e~^'^ 

CKM{2, 3) = [sin(a;„)Arf„ -f- cos(a;„)HrfJe"^‘*’ 

CKM{3,1) = [-cos{ud)A*^^ + sm{ud)BduW'^ 

CKM{3, 2) = [-sin(a;rf)A^„ - cos{ud)Bdu]e^'^ (53) 

where $ is the last of the hve that can be freely chosen. We chose a value for this so that 
the imaginary part of CKM{2,3) vanishes, corresponding most closely to the convention of 
the PDG|15| for the CKM matrix: 


<h = arctan 


sm(x)„)/m(/li„) 


(54) 



so that 


sin(<l>) 


sm{uJu)Im{Adu) 


(55) 


^J[cos{uu)Bdu + sm{uJu)Re{Adu)]‘^ + [sin(a;„)/m(Ad„)]2 
sm{uJu)Im{Adn) 


and 


cos($) 


yJcos‘^{uu)Bj^ + sin^ {uu)A*^^Adu + sm{2uu)BduRe{Adu) 

cos{uu)Bdu + sm{uu)Re{Adu) 


(56) 


Y/[cos(a;„)5d„ + sin(a;„)i?e(Ad«)]^ + [sin(a;„)/m(Arf„)]2 
cos{uu)Bdu + sm{uu)Re{Adu) 



which leaves 


CKM{2,3) = ^J[cos{uJu)Bdu + sin{uJu)Re{Adu)]^ + [sin(a;„)Jm(y4rf„)]2 (57) 


and accomplishes the desired resnlt. 

To leading order in e, CKM{3, 3) = 1, and the relation between uju and ojd is immediately 


hxed by the (2 x 2) Cabibbo rotation in the light qnark sector (the CKM{i,j) entries for 
i = 1,2 and j = 1, 2): 


COS^COd - UJu) Sm^LJd - UJu) 
-sin(a;rf - Uu) cos{ud - Uu) 

cos(0c) sin(0c’) 

-sin(0c) cos(0c) 


CKMc 


(58) 


(pins order e corrections), i.e., we identify 


(59) 


Qc = ojd - oju ■ 


At this point, it is apparent that the CKM mixing depends solely on the difference between 
the diagonalizations of the up-quarks and the down-quarks, as it must. 

A. PDG evaluation 

The PDG [l5] provides only the absolute values of the entries of the CKM matrix. It 
also presents a matrix form that has only real entries in the hrst row and third column of 
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the CKM matrix, except for the (1,3) matrix element. This is achieved by locating the 
one required phase in the matrix that produces rotation about the second axis, where the 
sequence of rotations is hrst about the third axis (almost Cabibbo), next about the second 
axis, and dually about the hrst axis, proceeding from right to left in the products in the 
usual way, viz. 



1 0 0 


cl3 0 sl3e-*^ 


cl2 sl2 0 

CKMpDG — 

0 c23 s23 

X 

0 1 0 

X 

-sl2 cl2 0 


0 -s23 c23 


-sl3e*^ 0 cl3 


0 0 1 


cl3cl2 cl3sl2 sl3e-*^ 

-cl3sl2 - cl2s23sl3e*>^ c23 cl2 - sl2 s23 sl3 s23cl3 

-cl2 sl3 - cl2 c23 sl3 -s23 cl2 - sl2 c23 sl3 c23cl3 


(60) 


where as usual, cl3 = cos( 6 ' 2 ), etc. and we have changed the PDG phase notation from 5 to 
X to avoid confusion with our mass ratio parameter above. 

We follow the PDG structure precisely below and our construction above agrees with its 
structure through hrst order in e as the sines of all of the angles are small (see below). 
However, we need to have explicit imaginary components for all of the matrix entries, rather 
than only moduli as reported by the PDG. To proceed, we have constructed a version of the 
PDG result where we assume that all three of the mixing angles reside in the hrst quadrant. 
This is not justihed, but demonstrates how the constraints on BSM parameters may be 
extracted were such information available. 

Taking values from the PDG [15], and using their parametrization, we obtain central values 
for the real and imaginary parts of these quantities in the relations: 


CKMl?> = 0.001067 - 0.003386/ 

CKM2?, = 0.04141 + 0.0/ 

CKM?,l = 0.008375 - 0.003294/ 

CKM32 = -0.04057 - 0.0007690/ (61) 

where the rhs in each case is the corresponding entry of the matrix CKMpdg when, as noted 
above, the particular set of signs for the sines is chosen corresponding to all three angles 
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being in the first quadrant. Also, using the entries in the upper left 2x2 block, (which are 
real through hrst order in small quantities as sl3 and s23 are both small) we estimate the 
value of the Cabibbo angle, ©c, as 


©c = 0.2291 


(62) 


i.e., approximately 13.2°. 

Other choices for extracting the full matrix elements could be investigated as well, but this 
demonstrates that at least one solution exists. We have investigated a number of alternatives 
and hnd that the largest differences, apart from signs, are in the real and imaginary parts 
of CKM13 and part of CKM32, but the changes are not large, e.g., ~ 20%. However, if 
the phase is placed in an alternate location, for example so that the hrst row entries are all 
real, then larger changes are obtained in the real and imaginary parts, although the moduli 
are maintained, of course. 

We uniformly present 4-digit values for consistency, but the changes between the 2012 and 
2014 reports suggest that in a number of cases the values are not known to better than two 
digits, at most, although some of the uncertainties are a small fraction of a percent. Due to 
the larger uncertainties, however, we conclude that carrying out our analysis to order is 
not warranted at this time. 


B. CP Violation 

We have taken advantage of the phase freedoms, particularly as noted by Kobayashi and 
Maskawa (KM) [T3], and used them to ensure that a large CP-violation, which would be 
inconsistent with experiment, does not arise within the (2 x 2) light quark sector. We must, 
however, examine what CP-violating implications are introduced by the BSM parameters 
that we have introduced that produce complex amplitudes. In particular, we can examine 
whether this is sufficient to be the only source of CP-violation. 

The invariant characterization of CP-violation was described by Jarlskog [19]. The Jarlskog 
invariant quantity, which we label J^ appears only at order e^, and is given by |T5| 

J = Xm[CKMijCKMk,iCKMliCKMlj] 
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( 63 ) 


= cos(6'i2)sin(6'i2)cos(6'23)sin(6'23)sin(6'i3)cos^(6'i3)sin(5) 

= (3.06 ±0.21) X 10“® 

up to an overall sign ambiguity, in the standard PDG representation of the CKM matrix. 

At the hrst order in e level of approximation, only the combination of matrix elements 
[i = j = 2, /c = Z = 3] reproduces the correct result for J\ 

J72233 = ±Xm[C'A:M(2,2)C'A:M(3,3)C'A:M*(2,3)C'A:M*(3,2)] 

= ±cos(0c) 1 Xm{[sm{uu)Adu ± cos(a;„)5rf„][sin(a;d)A^„ ± cos{ud)Bdu\} 

= ±cos{Qc)siia{Qc)BduXm{Adu) (64) 

We have checked that completing the unitary structure of Xtot to second order in e reproduces 
the correct result from any i, j, k, I combination. With the value of known, this provides 
one constraint on one pair of the combined parameters. 


C. BSM parameter constraints 


The parameter combinations in Eqs.( 53|[^52[59 ) display the fact that there are 4 quantities 
that can be directly related to the CKM matrix elements: Re{Aud)^ Irn{Aud)^ Bud and either 


one of Uu or Ud- Using the values of the CKM matrix elements in Eqs.(61), we can extract 
values for these combinations of BSM matrix elements. 


As seen above, ©c = Ud — oOu- If we further simplify the notation by dehning 


Bdu {^^dVcd ^uVcu) 

Y = Re{Adu) = edybdCos{Cd) - e«2/b«cos(C„) 

Z = Im{Adu) = edybdAn{Cd) - e„|/6„sin(Cn) (65) 


then in our leading approximation for the CKM matrix, the hrst two elements of the third 


column and row can be rewritten (making use of the phase freedom described at Eq.(53) 
and immediately following) as 


cos(2a;)XF - Isin(2n;) (X^ -Y^-Z^) iXZ 
CKM13 = - ^^^^ ± 


Q 


Q 


( 66 ) 
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CKM23 

CKM31 


Q (67) 

cos(c(;)sin(c(; + Qc)X'^ — sin(a;)cos(c(; + Qc){Y‘^ + — cos{2uj + Qc)XY 

Q 


+ 


iXZcos{ec) 

Q 


( 68 ) 


CKM32 = 


cos(a;)cos(a; + Qc)X^ + sin(a;)sin(a; + 0c)(F^ + Z^) + sin(2a; + Qc)XY 

Q 


iXZsm{Qc) 

Q 


(69) 


where we drop the index and write oj for oj^ and where 


Q = (cos(ci;)X + sin(c(;)y)^ + (sin(a;)Z)2 

= 0.04141 (70) 


and in these terms, 


J = ±cos(0c)sin(0c)XZ 


(71) 


It is straightforward to see from the imaginary parts that the (3,1) and (3, 2) elements 
contain no new information beyond that from the (1, 3) and (2, 3) elements, which is true 
for the real parts also. It is also clear that the imaginary parts are consistent with the form 


of J as given in Eq.(64). With a negative value for J we can solve for X, Y and Z as 
functions of Uu using the values of J, and of CKM13 and CKM23 in Eqs.(|M|. We hnd 


—0.00052189 + 0.0041169sin(a;„)cos(a;^l) — 0.070497cos^(a;„) 
•^2.9410 — 2.9194sin^(a;„) — 0.17049sin(ci;„)cos(ci;u) 
0.0020585 - 0.070497sin(a;„)cos(a;„) - 0.0041169cos2(a;„) 

1^2.9410 — 2.9194sin^(a;„) — 0.17049sin(a;„)cos(a;„) 
0.005729 

-y/ 2.9410 — 2.9194sin^(a;u) — 0.17049sin(a;„)cos(a;^i) 


(72) 

(73) 

(74) 


These functions are plotted in Fig.([^. Note that all quantities are no larger than of order 
ed (in fact, they are < 2ed) and so confirm that the BSM parameters are consistent with 
being “natural” in the usual sense. 


We note here that Eq.(70) demonstrates that Cartan sub-algebra perturbations alone are 
insufficient to satisfy these experimental constraints. If we set all of the j/j to zero except 
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for I/O, 2/3 and //s, then nsing Eqs.(2lp^40[44), we have 


1/3 ( 1 - 5 ) . ^ 

^ =ya= —^—sin(2n;) 

2/8 ( 1 - 5 ) , . 

71 — 2 — 


(75) 

(76) 


or 


2/3 ~ 

2/8 ~ 



sin(2a;) 

cos(2a;) 


(77) 

(78) 


Under these conditions, it also follows from Eq.(25) that yc = |/3^2/3. From Eq.(23), we 
similarly obtain yi, = |/8^2/3. Combining these in Q (with ^ = 0 so Z = 0 here), we hnd it 
bonnded by (nsing the larger e^) 


Qbound = < 0.02 (79) 

and so, far too small to match the valne in Eq.([70|) even if e„ were to contribnte positively. 
Note that nsing the phase freedom to prodnce CP-violation will not change this resnlt. 


Conversely, it is possible for the reqnired valne of Q to be attained with 1/3 = 0 or //g = 0, 


and perhaps even both. Even with Eqs.(80) below satished, this only reqnires 1/5 = 1/7 ~ 2, 
which is still “natnral”; 1/4^6 provides only a small contribntion. So, perhaps snrprisingly, 
not only are some non-Cartan BSM contribntions reqnired, it may be that both of the Cartan 
snb-algebra symmetry-breaking components (bnt not yo) may vanish and the entirety of the 
non-zero small masses and mixing may be dne strictly to non-Cartan BSM contribntions. 


D. Inverting Definitions 


It is possible to invert the seqnence of dehnitions of these qnantities and retnrn to the original 


yi in Eq.(19). The resnlt is somewhat cumbersome, bnt rednces the nnmber of independent 
parameters to be considered. Explicitly, we can set 


1/2 = 1/5 — yl and 
2/5 - 1/7 = 0 


( 80 ) 
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X. Y, Z vs. coup 



FIG. 1. Variation of combined BSM parameters for real and imaginary contributions to CKM 
matrix as functions of oju- 


These eliminate two of the imaginary terms in Eq. (19) so that this intermediate form of the 
mass matrix becomes : 


M 


intalt — 


+ e X 


0 0 0 
0 0 0 
0 0 1 

Ivo - \/ly8 + i(i/i - 2i/4p6) 

;yj(~2/3 + yime) 


(81) 


2/3 T yimo) \J^y^ y^p& ''/^ly'd) 

Jlyo + 7^2/8 - 2/1 (22/3 + 2/4m6) 

^32/8 T 372^^^^ 2/4p6 T V^f^2/9) yim &) 1 + 5 2^^220 + 3(2/1 + 2/4p6) 

where we have replaced the labels for the equal quantities 1/5 and yi by yg. (Recall that ^ 
does not appear as we hxed its value by the condition that the trace of Aiintait must be 
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equal to the trace of the (scale removed) diagonal matrix of eigenvalues.) 

Note that this procedure is not identical to the phase adjustments we made use of above. 
There the surviving real quantities in the mass matrix included contributions from the 
imaginary quantities through a normalization factor. In this truncated form, we can identify 


Va 

Vb 

Vc 




yime) 


3^2 
1 


1 /(2a/3|/8 + 2yi - + 36?/| 


ye 


(22/3 + yim .&) 


(82) 


and recall the unaltered yf = :^2/8 — |(22/i — 2 / 4 p 6 )- 


Using Eqs.(21, 40, 43, 44) and the assumptions in Eq.(80), we can simplify and invert those 
relations to obtain 


2|/i - yAp& = VSys “ ^(1 “ 5)cos(2cn) (83) 

\/3 

l/ 4 m 6 = 2/3 - ^(1 - 5)sin(2a;) (84) 

That is, we can rewrite these combinations of non-Cartan parameters in terms of j /3 and yg 
along with a dependence on oj. (As + ©c, this is still dependence on only a single 

angle, u = Uu-) Since only the combinations on the Ihs appear, the separate values of yi, 
j /4 and ye need not be specihed. This provides an example of how the parameters must be 
related in any BSM model. 


VI. DISCUSSION 

The most striking result of the view espoused here is that the smallness of the non-Cabibbo 
mixing follows from the ratio of the middle to largest masses of the quarks. Both are due 
to the perturbative size of BSM corrections to the initial “democratic” starting point. In 
contrast, the relatively large size of the Cabibbo mixing is allowed by the diagonalization 
process. A surprising result is that the BSM perturbations need not include all Cartan sub¬ 
algebra components, contrary to common analyses, while conversely, non-Cartan sub-algebra 
BSM perturbations are necessary. 
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The mass ratios of the quarks are scale dependent, and one could examine the effects of that 
scale dependence on the CKM matrix and our £t. However, even the ratios are generally 
not that well known and do not vary signihcantly with scale over the range from 2 GeV, 
where the lightest quark masses are generally dehned and determined, to the scale of the 
6-quark, nor from there to the weak scale which is also very close to the top quark mass. 
Rehnements responding to these issues are certainly warranted, but we do not expect them 
to produce large corrections to our BSM parameter constraints determined here. In fact, 
since the effects considered here are dominated by the values of e, only the uncertainties 
associated with the masses of the strange and charmed quarks should be significant, as 
the b- and t-quark masses are relatively accurately known. Fortunately, the very large 
uncertainties associated with the ratios of the two lightest quarks do not play a signihcant 
role in establishing the conhguration, although they will be important for precision analyses. 

We have carried out the straightforward extension of our results to the next higher order in e 
which might, in principle, be able to further constrain the values of the unknown parameters. 
Unfortunately, utilization requires knowledge of the relevant experimental values to order 
e^, i.e., to of order parts in 10"^, which is an accuracy generally not presently available. More 
accurate measurements would certainly change this conclusion. 


A. BSM contributions 

In general terms, Fig.([^ shows the nature of expected BSM corrections that could distin¬ 
guish the different fermions and lead to the small corrections that we hnd in our hts. The 
Lagrangian structure that we have in mind uses Weyl spinors for the separate left-chiral 
{dx, for member of a weak interaction doublet) and right-chiral {sx, for a weak interaction 
singlet) parts of the fermion Dirac bispinors, but nonetheless produces Dirac mass terms 
which may be simply represented as above. 

Interestingly, we expect these loop calculations to be hnite as they involve only differences 
within the triples of fermions. This effect was observed in Ref. [20] , where a symmetric, overall 
divergence appears but the differences in mass corrections are hnite, although that model 
is for a quite different application of symmetry-breaking mass corrections and employed 
vector gauge bosons. Although the calculation there was done only for Cartan sub-algebra 
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FIG. 2. A BSM loop correction to the fermion-Higgs-boson vertex that alters the charged 


mass matrix from “democratic” to that shown in Eq.(16). Left-chiral fermions with weak 
tions are labelled dx for “doublet” and right-chiral fermions with no weak interactions are 
Sx for “singlet” or “sterile”. 


fermion 

interac- 

labelled 


corrections, it is not unreasonable to expect that the off-diagonal corrections will also be 
hnite. Note that without the intermediation of the Higgs scalar vacuum expectation value, 
the dx is not changed to an Sx which prevents completion of the loop. 

Fig.([^ is drawn for a BSM scalar boson interaction, but a BSM vector could in principle 
also couple both the Weyl spinor dx to a Weyl spinor dx and similarly Sx to Sx- This simply 
requires interchanging the labels on either side of the Higgs’ coupling to complete the loop. 
Again, without the intermediation of the Higgs scalar vacuum expectation value, both the 
dx and the Sx pass through the loop unchanged without coupling to mass, so the BSM 
correction affects only vertex renormalization. 
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B. Leptons 


Of course, application to the leptons also comes to mind. For the charged leptons, Fig.([^ 
still applies, as it also does for Dirac mass terms for neutrinos. However, producing these for 
the neutrinos requires the existence of uncharged Weyl (right-chiral) fermions with no SM 
interactions at all, i.e., the so-called “sterile” neutrinos, as opposed to the known ones which 
are “active” with respect to the weak interactions. Under the long honored assumption 
of quark-lepton (Zweig-Glashow) symmetry, the existence of these states has been widely 
presumed since the early days of Grand Unihed Theories [ 21 ], starting with S'O(IO). Aside 
from the prediction of the charm quark, this symmetry (or regularity) successfully predicted 
the existence of the t- and 6-quarks as soon as the r-lepton was discovered. 122 ] There may 
even be some recent experimental evidence for “sterile” neutrinos. [23] 

Furthermore, with the discovery 121 of “Dark matter” in the Universe, it is clear that there 
are additional particles beyond those in the SM and that these particles are sterile. However 
many different types there are, only three distinct Weyl spinor combinations of these degrees 
of freedom are required to couple to the active neutrino Weyl spinors in conjunction with 
the Higgs boson to produce Dirac mass terms for the neutrinos. There is also no impediment 
for these sterile helds to acquire Majorana masses independently of the Higgs. 

It follows that a 6 X 6 structure in terms of Weyl spinors is required, rather than the simple 
3x3 Dirac mass structure that the charged fermions can be reduced to by the standard 
construction of Dirac bispinors. (Any additional sterile fermions can be block-diagonalized 
away, leaving the three needed, albeit perhaps of a complex structure in terms of the original 
model degrees of freedom.) The Dirac mass terms of the neutrinos, ttid, have the same form 
as that for the charged fermions as described above, but they now appear in 3 x 3 off-diagonal 
blocks in this 6x6 matrix. The upper left 3x3 block remains zero, as the SM does not 
produce Majorana masses for the active neutrinos, nor is it necessary for BSM physics to do 
so directly: The “see-saw” mechanism [6] produced by a sufficiently massive lower right 3x3 
block of “sterile” neutrinos leads to light, Majorana mass eigenstates that are dominated by 
active neutrino amplitudes. 

That lower right 3x3 block of Majorana masses (with a structure as given in Eq.([^ 
above) for the “sterile” Weyl spinors (corresponding to what would have been the right- 
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chiral component of a normal Dirac neutrino wavefunction) is unconstrained. As we pointed 
out many years ago |2S], neutrino mass mixing of the almost purely “active” eigenstates can 
be expected to be similar to that for the quarks, unless there is some particular structure to 
this 3x3 block of Majorana masses for the “sterile” Weyl spinors; i.e., barring any special 
circumstances, the leptonic analog of the CKM matrix, namely the PMNS-matiix in. 
should be similar to the CKM matrix. 

At that time, the concern was to determine whether or not neutrinos should be expected 
to have mass and whether or not their mixing should be expected to be large enough to 
measure. As we now know, the masses are very small but the mixing is even larger than 
for the quarks and very close to the particular TBM form that we showed above applies 
separately to the up- and down-type quarks, but cancels in the CKM matrix. 

We have identihed a6x6to3x3 block-diagonalization procedure that provides for a 
determination of the relations required between the sterile neutrino Majorana mass matrix 
and the BSM parameters in the neutrino Dirac mass matrix, so that, in the diagonalized 
6 x 6 mass matrix, there is negligible mixing of the active parts of the mostly active Majorana 
mass eigenstates relative to their initial structure. That is, for the leptonic weak interaction 
currents, we have ascertained that the conditions required, so that the factor contributed 
by neutrinos to the PMNS-maXiix is the identity (or close to it), may be satisfied. It 
follows that, under those conditions, the PMNS-matiix for the weak lepton currents will 
be almost, but not exactly, of the TBM form with small corrections certainly coming from 
the diagonalization of the charged leptons, which may even be the dominant corrections: 

PM NS ^ TBM. (85) 

This is consistent with current experiments which show that the PMNS-matrix is indeed 
quite close to the TBM-ma.tnx form, but the quantity ^ 13 , that vanishes in the exact limit, 
is not zero |26j . 

If the constraints determined in our 6 x 6 to 3 x 3 block-diagonalization procedure are 
satisfied, then the “democratic” plus BSM hypothesis for the fermion mass matrices would 
provide a unihed understanding of all of the weak current mixing structures simultaneously, 
subject to those additional constraints on BSM physics being satished. Conversely, any 


model of BSM physics that satishes these relations will produce the result in Eq.(85). We 



will present a detailed analysis of the leptonic sector in a future paper. 


VII. CONCLUSION 

We have started from the assumption that within the SM, for the fermions with a given 
electric charge, the Higgs doublet is sensitive only to the quantum numbers of the left-chiral 
Weyl spinor parts of the Dirac wave functions. With this assumption, a basis can be chosen 
so that the iso-singlet terms formed between them and the Higgs doublet couple equally to 
all of the right-chiral Weyl spinor parts. This in turn implies that the SM mass matrix should 
have a “democratic” form, with one massive eigenstate and two massless ones. Upon adding 
perturbative corrections of a completely general form, presumed to arise from BSM physics, 
we hnd that a consistent set of parameters may be extracted that conforms to the known 
quark mass spectra and CKM mixing matrix, including UP-violation. These parameter 
value constraints provide information on matrix elements of the manner in which BSM 
physics couples to SM degrees of freedom. The question of why the overall mass scale for 
the up-quarks is signihcantly larger than that for the down-quarks remains unresolved, but 
can be accommodated if there are a pair of Higgs bosons, as is required in supersymmetric 
models, for example. 

It is clear that, in this approach, extracting more detailed information on the nature of BSM 
physics and the value of BSM matrix elements requires a more accurate determination of 
the quark mass ratios and their mixing amplitudes in the weak interaction. It would be of 
great value if the separate real and imaginary parts of the CKM matrix elements could be 
determined experimentally. 

Under the “see-saw” assumption regarding the existence and nature of sterile neutrino com¬ 
ponents, the extension of these ideas to leptons can also produce the beginning of an under¬ 
standing as to both why the PMNS matrix is approximately of tri-bi-maximal form and also 
why it is not exactly so. Along with the information from the violations of “democracy”, 
the potential information on the sterile neutrino mass matrix may open the door to learning 
about the physics in the dark matter sector, with the sterile neutrinos as the hrst component 
known from other than gravitational interactions. 

Finally, we note that the intermediate propagation of sterile as well as active neutrinos in 
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variations of the graphs above applied to weak decay box graphs may be relevant to recent 
observations of violations of lepton universality, such as in Ref.f|27j). 
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